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The accurate analytical solution for the low temperature 1/f noise in a microwave dielectric 
constant of amorphous films containing tunneling two-level systems (TLSs) is derived within the 
standard tunneling model including the weak dipolar or elastic TLS-TLS interactions. The results 
are consistent with the recent experimental investigations of 1// noise in Josephson junction qubits 
including the power law increase of the noise amplitude with decreasing temperature at low temper¬ 
atures T < O.IK. The long time correlations needed for 1/f noise are provided by the logarithmic 
broadening of TLS absorption resonances with time due to their interaction with neighboring TLSs. 

The noise behavior at higher temperatures T > O.IK and its possible sensitivity to quasi-particle 
excitations are discussed. 
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1/f noise exists in a variety of physical systems 
and it dramatically restricts performance of modern elec¬ 
tronic and quantum nanodevicesi^^— The inverse fre¬ 
quency dependence of the noise power spectral density, 
^xxif) oc 1 //, is a consequence of a logarithmically slow 
relaxation. The slow dynamics is often associated with 
the random ensemble of fluctuators possessing a loga¬ 
rithmically uniform spectrum of relaxation timesi^i^i^^— 
Such fluctuators do exist in amorphous solids in the form 
of universal tunneling two-level systems (TLSs).— With 
the advent of superconducting quantum bits (qubits) 
based on Josephson j unctions a comprehensive study 

of the noise properties due to TLSs has become cru¬ 
cial for the achievement of high-fidelity quantum com¬ 
putation. TLSs are ubiquitous, appearing in wiring di¬ 
electrics, Josephson junction barriers and other disor¬ 
dered insulating regions. The deleterious effects of the 
coupling of TLSs to the qubit are associated with the 
resonant absorption of the qubit energy by TLSs in the 
microwave frequency range, as well as by 1/f noise in the 
microwave dielectric constant,— resulting in qubit deco¬ 
herence. 

Dissipation and 1/f noise in superconducting res¬ 
onators due to TLSs have been extensively studiedj^^— 
Recently, it was found that the 1/f noise amplitude in 
high quality superconducting resonators increases with 
decreasing temperature as with rj « O.3.— This 

dependence was considered as being incompatible with 
the standard tunneling model (STM)— 

A qualitative theoretical model interpreting this re¬ 
cent experiment has been proposed in Refs, [la and iisi . 
suggesting an energy-dependent TLSs density of states 
(DOS), g{E) Ri E^, in contrast with the standard TLS 
model suggesting a constant DOS— The authors used 
their assumption to interpret the additional exponent 77 
in the noise temperature dependence and the anomalous 
temperature dependence of TLSs decoherence rate dis¬ 


covered earlier using a single TLS spectroscopy— This 
assumption conflicts with earlier experimental data show¬ 
ing energy independence of the DOS, leading to logarith¬ 
mic temperature dependence of the dielectric constant 
and sound velocity2£i2i (see however Ref. . 

In this paper we derive an exact analytical solution for 
the low temperature 1/f noise in a microwave dielectric 
constant of amorphous films, assuming a homogeneous 
DOS as in the STM. This solution can be directly com¬ 
pared with the available experimental data and used to 
quantitatively characterize TLSs properties which is im¬ 
portant for understanding of the nature of TLSs and the 
reduction of their destructive effects in superconducting 
quantum devices. Particularly, we show that by properly 
considering the logarithmic temperature dependence of 
the noise power spectral density, the dependence 

of the amplitude of 1/f noise on temperature can be 
explained within the assumptions of the STM and pre¬ 
dict the noise dependence on the external field amplitude. 
The derivation focuses on the relevant regime of low tem¬ 
peratures, where the thermal energy is much smaller than 
the field quantization energy, namely k^T huj. 

At higher temperature, k^T > hu, we argue that TLSs 
based models cannot explain the noise temperature de¬ 
pendence. An alternative noise source in this regime 
is discussed by employing the possible contribution of 
quasi-particles excitations in the superconductors. It is 
demonstrated that such excitations may be responsible 
for the anomalous temperature dependence of the relax¬ 
ation and decoherence rates of TLSs, observed in this 
regime— 

Consider the TLSs contribution to the dielectric con¬ 
stant measured using the input field, applied along 

the z axis at a microwave frequency to and low tempera¬ 
ture kgT ^ hu. In this regime the resonant contribution 
of TLSs to the real part of the dielectric constant domi- 
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nates and can be expressed in the 
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herence rate is composed of the contributions of relax¬ 
ation and pure phase decoherence, T^Y = (2T’ii)~^+T'“^, 
where the phase decoherence rate T^l is determined 
by the TLS spectral diffusion induced by its interac¬ 
tion with neighboring thermal TLSs, i.e. TLSs for which 
E ^ Aq m This rate is given b y^^i^^ 


Here, V and e represent the sample volume and bare 
dielectric constant. The summation is taken over all 
TLSs i having transition dipole moments Pj^j, relax¬ 
ation and decoherence times and T 2 j, respectively, 
and time-dependent energies E^t). The z component 
of the transition dipole moment is related to the corre¬ 
sponding component of the TLS dipole moment, p^^^, by 
Ptr,iz = Piz\t/Fir^ where Ap, is the tunneling am¬ 
plitude of TLS i. Following the STM^i^ we assume 
that TLSs possess the universal distribution P{Aq,E) = 
PqE/{Aq^JE'^ — Aq) with respect to their energies and 
tunneling amplitudes, where Pg is a material dependent 
constant. 

The time dependence of the energies is induced by the 
spectral diffusion caused by weak dipolar or elastic TLS- 
TLS interactionj^^ The time dependence of the energy 
can be treated classically because at the temperatures 
under consideration, k^T <C huj ~ E, the equilibrium 
thermal fluctuations can induce transitions of TLSs only 
with an exponentially small probability 

Relaxation of TLSs is caused by phonon emission or 
absorption, with the relaxation rat o^°’^^i^^ 
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where the proportionality constant A ~ 10®s ^ char¬ 

acterizes the TLS-phonon interaction.— The TLS deco- 


1 YMxkBT-AT^ 

Ti- = V X-S- ’ 


where A^ = ^Ef — Ag^ is the TLS asymmetry and the 
dimensionless constant x = ^ ’ 10~^ represents 

the universal product of the TLSs DOS and their 1/R^ 
interaction strength.— The product X^bF represents the 
typical interaction with thermal TLSs and AT^ repre¬ 
sents the relaxation rate of thermal TLSs. Although the 
use of Eq. ([3]) in the expression for the dielectric constant 
[Eq. H])] has not been justified theoretically, its approxi¬ 
mate relevance was demonstrated experimentally— One 
should notice that in the linear regime of weak external 
field the related term is negligible in the noise spectral 
density (see Eq. ([5]) below and estimates of parameters 
of this equation in the next two paragraphs). In the op¬ 
posite regime of a very strong held the accuracy of our 
results is limited to the above assumption and may need 
a special consideration. 

The noise is determined by the time-dependent cor¬ 
relation function of the dielectric constant huctuations 
SxxY) = ('^eTLsW'5eTLs(0))/e^ where SctlsY) = 
f-TLsit)~ Y-tls)■ Correlations between different resonant 
TLSs contributing to the noise can be neglected because 
of the weakness of the TLS-TLS interactions (see e.g., 
Ref. [23). The correlation function can then be expressed 
as 
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The correlation function in Eq. o can be evaluated 
by means of the distribution function D{E\EqA) for the 
TLS energy E at time t, assuming E{0) = Eq at t = 
0— For 1/i?^ interaction, this distribution function is a 
Lorentzian 


D{E\Eo,t) 


1 W{t) 

Yw^{t) + {E- EqY’ 
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with a characteristic width W{t) = IFg(t)|A|/E, where 


VFg(t) 
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and T^^{y) = 8AT^y^ cothy. 

Low frequency 1 /f noise is determined by long times 
t > Is ^ (AT^ 1 ~^ ll^ where {ATY~^ estimates the min¬ 
imum relaxation time of thermal TLS. In this limit the 
width of the distribution Wg(t) grows logarithmically 
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with time and can be approximated as 

(7) 

This logarithmic increase of the resonance width is re¬ 
sponsible both for the appearance of 1// noise and for 
the additional exponent r] in the temperature dependence 
of the noise amplitude. This dependence has not 
been considered in Ref. [H, leading the authors to the 


conclusion that the experimental results are inconsistent 
with the STM. 

At low temperatures W{t) <C hw and one can per¬ 
form averaging of the correlation function in Eq. over 
the distribution of energy fluctuations [Eq. ([S])] and over 
initial TLS energies within the resonant approximation. 
Straightforward evaluation of integrals with respect to 
E{t) and E{0) yields 


.%Jt) = tanh^ 


47r3Po dx I 

hVe^ Jo x{l - x^) Jo 


where x = Ao/huj and y = cosO. Averaging in Eq. 
dS]) is made over absolute values of TLS dipole mo¬ 
ment, Po, forming an angle 6 with the AC field. In 
the numerical calculations below we assume Po ~ 5D, 
which is well justified experimentallyii The relaxation 
rate T^^{x) is given by Eq. m with Ej replaced by 
huj and Apj/Ej replaced by x. Similarly, T^^ix) = 
{2T-i^{x))~^ + T~^{x), where T~'^{x) is given by Eq. ([3]) 
with IA^|E^ = y/l — (Aoi/Ei)^ replaced by \/\ — x^. 

The power spectral density of noise, can be 

evaluated as a Eourier transform of [Eq. ([8])] in 

the low frequency limit fT^, f /W ^ 1. It has the pure 
1// spectrum if the function S^^(t) depends on time as 
A — B\n\t\, which has a Fourier transform 5/(2/) at 
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/ ^ 0. The correlation function S^^(t) can be repre¬ 
sented using the expansion near tj Ri l/(27r/) in the 

approximate form 

(higher order expansion terms are smaller by the factor 
ln~^(Ar^//) r; O.I for the small frequency of interest 
/ R! 0.1 Hz). The 1/f noise power spectral density can 
then be expressed as = —{l/2f)dS^^{tj)/d\ntf. 

For a quantitative comparison of the theory with the 
experiment it is convenient to introduce a volume inde¬ 
pendent parameter in a similar way to the experimentally 
determined Hooge’s constant for 1/f conductivity noise 
in semiconductors)^ We define this parameter as 
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which is the ratio of the noise spectral density multiplied 
by the number of thermal TLSs, Nrp = PoVk^T, to the 
squared average loss tangent due to the TLSs, tani5 = 
(e")/e = (47r^/3e) PQ(pg) tanh (fiw/2fcBr). One should 
notice the logarithmic weak frequency dependence of the 
noise amplitude equivalent to a hardly distinguishable 
dependence for the typical system parameters (see 
the caption to Fig.[T|). 

Consider the noise amplitude temperature dependence 
in the linear regime > 0. For typical TLS pa¬ 
rameters X ~ 5 • 10“'^ and A r; one has 

W'o(r) - 10®r(s-i), r; 10®T2(s- 1) and Tf^ r; 

10®s“^. Consequently, the second term in the denomi¬ 
nator of the integrand in Eq. m can be approximately 
neglected at temperatures 0.01K< T < O.IK. The tem- 


3Wg(l/(2^/)) -f VI + {xypoF^cyT,{x)T^{x)/h^ 


perature dependence of the noise amplitude can thus 
be approximated by S^^{f) oc T~^ ln“^ (3.3 • AT^/27r/). 
The approximate exponent in the power law tempera¬ 
ture dependence oc observed in the exper¬ 

iment can be estimated as ?7 = — 1 — (ilnS' 2 , 2 ,(/)/(ilnT = 
—dhiarpj^g/dhiT = 6/ In (3.3 • AT^/27r/).— Setting / = 
O.lHz^ we obtain rj Ri 0.42, in agreement with the 
experimental observation. The actual dependence is 
slightly weaker because of the logarithmic integral in Eq. 
(|ni). Thus, the reported anomalous temperature depen¬ 
dence of the 1/f noise amplitude in superconducting res¬ 
onators at low temperatures can be interpreted within 
the framework of the standard tunneling model. 

The temperature dependence of the noise parameter 
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FIG. 1: The exponent 77 in the noise temperature depen¬ 
dence a-rj.s oc T”'' and the noise parameter (inset) 

vs. temperature. TLS parameters are chosen as^ po = 5D, 
PoUrms = 5 • 10-^ A = 10®s-lK-^ 10^2%- 5GHz. 


Otls ( 0 )] and the exponent 77 extracted from it are 
shown in Fig. [T] for various external field amplitudes. It 
is evident that in the linear regime, Fac 0 , the tem¬ 
perature dependence observed in the experiment in the 
regime k^T tuo is qualitatively reproduced. The ac¬ 
curate comparison of the theory with the experiment is 
the subject of a separate work that will help to determine 
the system parameters from the optimum data fit. 

It should be emphasized that our results do not rule out 
the possibility of some energy dependence of the DOS. 
For the specific case of 1// noise it shows that a signifi¬ 
cant contribution to the exponent 77 may arise due to the 
broad spectrum of relaxation times in amorphous sys¬ 
tems, even if the DOS is energy independent. However, 
other deviations from the STM, such as the non-integer 
exponents of the specific heat and thermal conductivity 
may originate from an energy dependence of the DOS at 
low energies. Further investigation is thus necessary in 
order to shed light on the energy dependence of the DOS. 
This can be achieved by additional measurements, such 
as temperature and frequency dependent dielectric losses 
or internal friction in the plateau regime, which should 
be both sensitive to the energy dependence of the DOS 3 

The strong reduction of the noise with increasing held 
amplitude, as well as the weakening of its temperature 
dependence are clearly seen in Fig. [1] These hndings 
are in agreement with the experimental resultsYet, 
in the large held limit, F^q > {hWQ{t)/pQ) ~ 

lOOOV/m, Eq. (0) predicts oc TF^q, while the ob¬ 
served held dependence is weaker This might be due to 
nonequilibrium heating^ and insufficiently large AC held 
used experimentally. A detailed analysis requires under¬ 
standing of the mechanisms of heat exchange between 


the superconducting circuit and the environment, which 
is beyond the scope of this letter. 

At high temperatures, k^T > hw, Eq. m as well 
as the earlier wor k^^d^ predict the strong reduction of 
the noise with increasing temperature due to the factor 
tanh^ {E^/2kBT) oc T~^. Such temperature dependence 
does not ht the experimental observation^ which re¬ 
ports a weaker temperature dependence. We expect that 
different mechanism may be responsible for the 1// noise 
in this regime. According to our preliminary analysis this 
cannot be the relaxational interaction of TLS with the 
AC field^S since it does not lead to 1// noise. Indeed, 
the low frequency noise should be due to slow TLSs with 
the large relaxation times ~ l/27r/. They can in¬ 
deed contribute to the fluctuations of a zero frequency 
dielectric constantiS while their response to an AC held 
vanishes as f/uj^ 

Alternatively, the noise in dielectric tunneling bar¬ 
rier incorporated in superconducting devices may be 
due to quasi-particle excitations. In particular, a re¬ 
cent study has revealed an anomalous temperature de¬ 
pendence of the relaxation and decoherence rates of sin¬ 
gle TLSs in the amorphous tunnel barrier of a super¬ 
conducting phase qubit, showing a deviation from the 
single-phonon relaxation rate [Eq. (0)] at temperatures 
above O.IKJ^ As shown in Fig. 0 the relaxation rate 
of these TLSs can be Htted by combining the contribu¬ 
tions to the relaxation due to TLS-phonon interaction 
and the thermally activated quasi-particle contribution^^ 
as Tj“^(r) = Acoth (E/2knT) + . For both 

TLSs studied in Ref. [l^ the htting parameter A is close 
to the A1 superconducting gap 3.95fc^K. These quasi¬ 
particles may also be responsible for 1// noise at tem¬ 
peratures exceeding O.IK. 



FIG. 2: Experimental datai^ for temperature dependent re¬ 
laxation times of two different TLSs interpreted assuming the 
presence of quasi-particles. 

In conclusion, we calculated the I// noise in mi¬ 
crowave dielectric constant produced by TLSs in amor¬ 
phous solids. The STM involving the long range TLS- 
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TLS interactions has been used. The results are consis¬ 
tent with the recent experimental data at low temper¬ 
ature, whereas at higher temperature other mechanism 
may be responsible for the noise, possibly one associated 


with broken Cooper pairs. 
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